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For those angular multipoles where cosmic variance is an issue, non-Gaussianities in the Cos- 
mic Microwave Background (CMB) anisotropics will be hard to detect. Here, we construct 
explicitly the best unbiased estimator for the CMB angular bispectrum. 
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CO . While the statistical properties of the CMB anisotropics are a powerful means to discriminate 

amongst the possible cosmological scenarios, actually measuring non-Gaussianity in the data is 
a very difficult taslfi. The typically small signal should be compared to the noise and the key 
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quantity is the signal to noise ratio. The noise creeps into the dataset through instrumental 
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errors, foregrounds contamination or incomplete sky coverage. Add to this the so-called 'cosmic 



variance': the fact that we only have access to one realization of the temperature anisotropics 
A( e ) = 5T/T(e) whereas theoretical predictions are expressed through ensemble averages. It 
can dominate the other sources of error and therefore, if one wants to unveil non-Gaussianity, it 
is necessary to address the cosmic variance problem for those quantities characterizing a possible 
non-Gaussian CMB temperature anisotropy distribution. For that one constructs estimators by 
performing spatial averages on the celestial sphere and finds the one which has the smallest 
IT/ \ possible variance. We here show the best unbiased estimator for the angular bispectrum Q 1 ^ 2 £ 3 

and we display the corresponding cosmic variance as well. Recall that Cg x g 2 t> 3 contains all the 
information available in the three-point correlation function, or its variants, like the skewness, 
collapsed and equilateral configurations! The present analysis borrows from recent joint work 
with Jerome Martml'i, to whom I am greatly reconnaissant. 
Expanding the anisotropics over the microwave sky as usual 

A(e) = J2a?Y f m (e) (1) 

im 

the first three moments can be written as 

(af) = 0, {a^a^*)=C h 5 hl2 5„ 2 , <a£»™ 3 ) = C ^ ^ 

where is a Wigner 3j-symbol. The second equation ensures the isotropy of the CMB. 

The quantity (a™ 1 ^ 2 *) is the second moment of the a™'s and Ci is usually called the angular 
spectrum. In the third equation, the quantity (a™ 1 a 7 ^ 2 alj^ 3 ) is the third moment while C^ 2 £ 3 is 
called the angular bispectrum. The presence of the 3j-symbol guarantees that the third moment 
vanishes unless mi + m<i + 1713 = and \li — £j\ < it < &i + £j- Moreover, invariance under 
spatial inversions of the three-point function implies the additional rule Em i\ + £2 + ^3 = even, 



in order for the third moment not to vanish. Finally, from this last relation and using standard 
properties of the 3j-symbols, it follows that the angular bispectrum is left unchanged under any 
arbitrary permutation of the indices £{. 

Let us call £ (C444) the estimator for the angular bispectrum £4^4- The most general 
definition reads 

£{C tMa ) = J J Jdn 1 dn 2 dn 3 E e /^(e 1 ,e 2 ,e 3 )A(e 1 )A(e 2 )A(e 3 ). (3) 

where E^ 1 ^ 3 is the weight function. The angular bispectrum is a real quantity and so is 
its estimator. Therefore, the weight function can be taken real. It is also symmetric under 
arbitrary permutations of directions ej. In addition, like Ce 1 e 2 e 3 i the weight function satisfies 
Eg 2 3 = E^ 1 ^ 3 , as well as for any other arbitrary permutation of the indices £{. The weight 
function can be expressed on the basis of the spherical harmonics as 

E^Hene^s) = E E E df/VV f^^Y^Y^). (4) 

i'^m!^ ^2 m 2 ^3 m 3 I m 'l m 2 m 3 / 

The properties of the weight function imply that the coefficients d must satisfy 

^ h (2 i 3 * _ ^i^m^+m'z+m^ «i <2 l 3 ^ «i «2 I3 _ ^ 'i <?2 H 

J *1 ^2 ^3 L J ^1 ^2 ^3 L J ^1 ^2 ^3 I J £ 2 ^1 ^3 l 

I m 1 j | — ~ m 2 ~ m 3 f I m l m 2 m 3 I 1 m 2 m l m 3 J 

where the last relation is in fact valid for arbitrary permutations of any two columns of the 
collective subindex. Like the weight function, d is also left invariant under arbitrary permutations 
of indices £i (not primed). The estimator can be expressed in terms of the coefficients d and 
the a£P's only: inserting the expansion of the weight function in the above expression for the 
estimator and using standard properties of the spherical harmonics one obtains 

Wm) = E E E df\ \ \ x <H?< 3 - (6) 

i^m^ ^2 m 2 ^3 m 3 1 m l m 2 m 3 / 

In practice, CMB observational settings are devised such that both the monopole and the dipole 
are subtracted from the anisotropy maps. This means that the coefficients d in the last equation 
are only non-vanishing for indices £[ > 2 in the collective subindex. Moreover, the coefficients 
d satisfy ^1+^2 + ^3 = even. We must now require that our general estimator given by Eq. 
@ be unbiased, i.e. (£(C&xi 2 h)) = Cfifefe- This forces the coefficients d to fulfill the following 
constraint 
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where we have defined a symmetrized Kronecker symbol for the £ multipole indices only, as 
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follows 5 g l 3 = \i3i 1 i l ^>i 2 i'J>i ?> i> i + 5 additional permutations ). It is easy to check that the con- 
straint equation satisfies the conditions imposed by Eqns. (||) on the coefficients d. Using the 
previous properties for d, relabelling the indices m'i <-> m' 2 in Eq. (^) and finally noting that 
l'\ + £ 2 + £'3 = £\ + £ 2 + £3 = even, which allows us to permute any two columns of the Wigner 
3j-symbol, one verifies that the left hand side of the constraint is invariant under £[ <-> £' 2 . The 
same applies for any pair of £ multipole indices and this explains the presence of the symmetrized 

6g 3 in Eq. ([?]) . We see from this that all coefficients d that do not satisfy £' x + £' 2 + £'3 = even do 
not enter the constraint. These terms only increase the variance (which we want to minimize) 
and as a consequence one can take them equal to zero. 

We are now in a position to calculate the variance of the estimator. Looking at Eq. (||) we 
see that this requires the computation of the sixth moment of the a™'s. After having made use 



of the properties of the coefficients d and rearranging the resulting 15 terms into two groups, 
straightforward algebra yields 



([£( c *l*2* 3 )] 2 ) - E E E fyCqpt 

e' 1 rn' 1 l' 2 m! 2 e' 3 m' 3 



(8) 



6dt\/\/ 3 * d 1 , 1 ,*,' 3 , s +9(-iri- 

I m-y m'2 m 3 J I m 2 J 



J 1 1 ^3 L J ^2 ^2 ^3 I 

J — m.g J | m 2 — ro 2 J 



The square of the variance of £(£^4) is given by °f(Q iJfa< ) = ([f (C^^)] 2 ) - (£ (C^^ 3 )) . 
Since departures from Gaussianity are expected to be small (specially on large angular scales), 
higher moments will be calculated in the mildly non-Gaussian approximation. Within this 
approximation we can write a™ = a^ ^ + ea™' 1 ' + 0(e 2 ) where a™^ is a Gaussian random 
variable and the expansion parameter e is small. [The label will be dropped out hereafter] 
Now, the term ([£(Cp 1 i 2 p 3 )] 2 ) is of order e° whereas the lowest non-vanishing order of (£(Cp 1 p 2 p 3 )) 2 
is e 2 . Therefore, the latter one will not enter the minimization procedure and the variance 
squared will be written as a% c « ([<S(Q 1 £ 2 £ 3 )] 2 ). This does not occur for the two-point 

correlatoEl; in that case both terms contributing to the square of the variance are of the same 
order in e . 

With a bit of effort one can see that the various contributions of the imaginary part of the 
coefficients d to the two terms, 6d*d and 9d*d, only increase the variance. Since we know that a 
vanishing imaginary part does satisfy the constraint Eq. (0), it can be disregarded in the sequel. 
Therefore, Eq. (0) can then be written solely in terms of real coefficients d as follows 
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Our next move now is to minimize this variance with respect to the coefficients d, taking into 
account the constraint of Eq. (0) 



+ A ii44 d r 4 4 4 1 1 m ' 



0. (9) 



Performing the variation 5 having in mind that the symmetries of the coefficients d must be 
respected, we get 
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" 1' -> 2' " 
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[. . .] terms are shorthand for the first one on the second line. This formula, together with Eq. 
(0), form a set of equations which completely determines the best unbiased estimator. 

From this last equation and using the constraint Eq. (^) we can get the general expression 

for the Lagrange multipliers. Thus, we multiply Eq. (|i~0| ) by the appropriate 3j-symbol and we 
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sum over the three indices m i . The first term is exactly the constraint and produces a o s 3 . 
Using the fact that a triple sum over the m^'s of the squared of a 3j-symbol gives unity, the 
second term yields the Lagrange multipliers themselves. Unfortunately, I don't have enough 
space to show that the last three terms vanish. Then, the Lagrange multipliers are given by 
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Plugging this into Eq. flTo|), one has 
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(12) 
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This is the final equation to be solved in order to determine the best unbiased estimator. A 
solution is 

'2 '3 _ / t'l 
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which leads to 



£Best(C£i^3/ 




m'i m'2 m'n 



1 1 2 3 

a h a h a h ■ 



(13) 



(14) 



Seems familiar? An estimator restricted to the diagonal case l\ 
to £2 = ^1 + 2 and £3 = i\ — 2) has been proposedQ'B for Bi = C/. 
was not to seek the best estimator, but to use the correspo 
non-Gaussian features of the 4-yr COBE-DMR data (see alscB 



-2 = £3 (and then extended 
j] The aim of these authors 
say, £(Bi) to analyse the 
1 L 1 1 * While their estimator is 



not unbiased, for it does not satisfy the constraint (^), by just removing an overall prefactor one 
gets our best unbiased estimator £Best(C^ 2 ^ 3 ), Eq- (0)- 

We now know the best unbiased estimator for C^i^ and then we can compute its variance, 
the smallest one amongst all possible estimator variances, which yields 



>) 



C^C^C^ (1 + 5^ + 5e 2 e 3 + 5^ + 2 <5^ 2 <5. 



(15) 



We like to dub this (the square of) the 'bispectrum cosmic variance' in perfect analogy with 
2Cj/(2£ + 1), which is (the square of) the variance of the best unbiased estimator 



for the angular spectrum, commonly known as the 'cosmic variance'. 
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a B like Bispectrum (like speCtrum) . . . or C, B, Al, .... Hope the Trispectrum will be called T. 
Still, since no ambiguity arises I stick to Ce 1 e 2 e 3 (for now). 



